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Abstract 
A  complex  matrix  is  said  to  be  stable  if  all  its  eigenvalues  have  negative 
real  part.  Let  7  be  a  Jordan  block  with  zeros  on  the  diagonal  and  ones  on  the 
superdiagonal,  and  consider  analytic  matrix  perturbations  of  the  form  A(c)  = 
J  +  cB  +  0{(^),  where  c  is  real  and  positive.  A  necessary  condition  on  B  for 
the  stability  of  A(()  on  an  interval  (0,fo),  and  a  sufficient  condition  on  B  for  the 
existence  of  such  a  family  A{e),  is  (i)  Re  tr  5  <  0;  (ii)  the  sum  of  the  elements 
on  the  first  subdiagonad  of  B  has  nonpositive  real  part  and  zero  imaginary  part; 
(iii)  the  sum  of  the  elements  on  each  of  the  other  subdiagonals  of  B  is  zero.  This 
result  is  extended  to  matrices  with  any  number  of  nonderogatory  eigenvalues  on 
the  imaginary  axis,  and  to  a  stability  definition  based  on  the  spectral  radius.  A 
generalized  necessary  condition,  though  not  a  sufficient  condition,  applies  to  arbi- 
trary Jordan  structures.  The  proof  of  our  results  uses  two  important  techniques: 
the  Puiseux-Newton  diagram  and  the  normal  form  of  Arnold.  In  the  nonderoga- 
tory case  our  main  results  were  obtained  by  Levantovskii  in  1980  using  a  different 
proof.   Practical  implications  are  discussed. 

1      Introduction 

A  matrix  is  said  to  be  stable  if  all  its  eigenvalues  lie  in  the  open  left  half  of  the  complex 
plane.  A  problem  of  great  theoretical  interest  and  practical  importance  is  the  following: 
given  a  matrix  A  on  the  boundary  of  the  set  of  stable  matrices,  i.e.  with  no  eigenvalue 
having  positive  real  part  and  one  or  more  eigenvalues  lying  on  the  imaginary  axis,  what 
perturbations  to  the  matrix  are  associated  with  stability?  To  be  more  specific,  suppose 
that  yl'"'  is  an  n  X  n  Jordan  block  Jo,  with  zeros  on  the  diagonal  and  ones  on  the 
main  superdiagonal,  and  consider  a  real  perturbation  A(€),  with  A{0)  =  Jq.  It  is  well 
known,  e.g.  [25,  Section  2.2],  that  if 


AU)  = 


0     1 


then  the  multiple  zero  eigenvalue  splits  into  n  eigenvalues 

i^l^u",    h=  l,...,n, 

where  w  is  the  nth  principal  root  of  unity  exp(27r!/n).  Thus  if  n  >  2,  at  least  one  of 
these  eigenvalues  has  positive  real  part  for  all  f  /  0. 

In  this  paper  we  give  a  much  stronger  result.  We  show,  for  example,  that  if /1(f)  is 
a  real  analytic  matrix  function  of  the  form 

^(e)  =  Jo  +  fB  +  0(e-) 

then  a  necessary  condition  on  B  for  the  stability  of  A{i)  on  an  interval  (0,fo),  and  a 
sufficient  condition  on  B  for  the  existence  of  such  a  family  A{c),  is 

(i)  tr  5  <  0,  i.e.  the  sum  of  the  diagonal  elements  of  B  is  less  than  or  equal  to  zero; 

(ii)  the  sum  of  the  elements  on  the  first  subdiagonal  of  B  is  less  than  or  equal  to  zero; 

(iii)  the  sum  of  the  elements  on  each  of  the  other  subdiagonals  of  B  is  exactly  zero. 

Note  that  (iii)  excludes  (1)  for  n  >  2.  In  fact,  we  show  that  the  jih  diagonal  in  the 
lower  triangle  (counting  the  main  diagonal  a.s  the  first)  is  associated  with  perturbations 
of  the  form 

where  w  is  the  jth  principal  root  of  unity.  This  largely  explains  (i)-(iii),  since  eigenvalue 
splitting  separated  by  angles  of  less  than  tt  in  the  complex  plane  cannot  be  tolerated 
for  stability,  while  the  splitting  associated  with  j  =  2,  namely  in  opposite  directions 
in  the  complex  plane,  is  permitted  only  along  the  imaginary  axis.  The  case  j  =  \ 
corresponds  to  an  eigenvalue  shift,  permitted  only  in  the  negative  direction.  Our  results 
are  established  in  quite  full  generality,  allowing  yl*°'  and  A(e)  to  be  real  or  complex 
and  requiring  knowledge  only  of  the  Jordan  form  of  /l'°'.  In  the  case  of  a  complex 
perturbation  of  Jo  (taking  e  to  be  real),  we  obtain,  instead  of  (i),  that  Re  tr  B  <  0, 
and,  instead  of  (ii),  that  the  sum  of  the  elements  on  the  first  subdiagonal  has  real 
part  less  than  or  equal  to  zero  and  imaginary  part  equal  to  zero.  This  is  because  an 
imaginary  shift  in  the  eigenvalues  is  permissible,  but  even  a  small  rotation  of  a  split 
along  the  imaginary  axis  is  not. 

In  the  case  that  j4'°'  has  one  or  more  derogatory  multiple  eigenvalues,  we  obtain 
a  generalized  necessary  condition  but  no  sufficiency  condition.  For  example,  in  the 
nondefective  (semisimple)  case  (e.g.  /l'°'  =  0),  conditions  of  the  form  (ii)  and  (iii) 
vanish  (since  a  perturbation  of  a  nondefective  eigenvalue  is  always  Lipschitz),  and  (i) 
remains  an  obvious  necessary  but  insufficient  condition  for  stability.  In  fact,  a  stronger 
necessary  and  sufficient  condition  is  well  known  for  the  nondefective  case;  if  j4'°'  =  0, 
it  is  trivially  that  the  maximum  real  part  of  the  eigenvalues  of  B  be  less  than  or  equal 
to  zero.  For  a  more  complete  discussion  of  the  semisimple  case,  see  [14,22].  To  a  large 
extent  the  present  work  arose  through  attempting  to  answer  the  questions  raised  in 
Section  6  of  [22].  The  present  paper,  together  with  related  work  in  progress  [4,5],  not 


only  answers  those  questions,  but  establishes  much  stronger  results  than  anticipated 
by  [22].  For  the  much  more  special  symmetric  or  Hermitian  case,  see  [21]. 

The  proof  of  our  results  uses  two  powerful  techniques:  the  Puiseux-Newton  diagram 
for  describing  roots  of  polynomials  with  analytic  coefficients,  and  the  normal  form  of 
Arnold  for  perturbation  of  multiple  eigenvalues.  These  are  both  explained  in  the  text. 

In  the  nonderogatory  case  our  main  results  have  already  been  obtained  by  Levan- 
tovskii  [16,17,18],  who  used  a  different  technique  to  study  C'  (continuously  differen- 
tiable)  perturbations.  Levantovskii  also  made  use  of  the  Arnold  normal  form,  but  he 
did  not  use  the  Puiseux-Newton  diagram,  which  we  consider  to  be  a  key  point  for  the 
understanding  of  these  results.  (Analyticity  is  required  for  the  use  of  the  diagram.)  See 
Levantovskii's  papers,  as  well  as  [2,  p. 255-257],  for  some  further  results,  particularly 
illustrations  of  the  shape  of  the  stability  region  in  parameter  space  for  various  cases. 

We  note  that  the  Puiseux-Newton  diagram  has  been  used  in  a  related  context, 
namely  the  stability  of  methods  for  solving  ordinary  differential  equations;  see  [12]. 
However,  that  work  has  a  very  different  emphasis,  not  being  concerned  with  multiple 
eigenvalues  but  instead  with  second-order  effects  of  the  perturbation  of  simple  eigen- 
values. 

All  the  results  have  a  straightforward  extension  to  the  case  of  reducing  the  spectral 
radius  of  a  matrix  with  multiple  eigenvalues.  This  case  is  relevant,  for  example,  when 
the  definition  of  matrix  stability  requires  the  magnitude  of  the  eigenvalues  to  be  less 
than  one.  This  situation  typically  arises  in  practice  when  studying  the  solution  of 
difference  equations  rather  than  differential  equations. 

The  paper  is  organized  as  follows.  Section  2  uses  the  Puiseux-Newton  diagram  to 
study  the  roots  of  polynomials  with  complex  coefficients.  Section  3  uses  the  Arnold 
normal  form  to  extend  these  results  to  perturbation  of  a  complex  matrix  with  multiple 
eigenvalues,  both  the  nonderogatory  and  derogatory  cases.  Section  4  briefly  discusses 
the  real  case.  Section  5  is  concerned  with  reducing  the  spectral  radius  of  a  matrix. 
Section  6  discusses  practical  implications. 

2      Roots  of  Polynomials  with  Complex  Coefficients 

A  complex  function  0(c)  is  said  to  be  analytic,  equivalently  holomorphic,  near  f  =  0 
if  it  may  be  expanded  in  a  power  series  in  e,  convergent  in  a  neighborhood  of  f  =  0. 
Consider 

P(A,f)  =  (A  -  Ao)"  +  /?i(0(A  -  Ao)"-'  +  . . .  +  /?„(0  =  0,  (1) 

a  polynomial  equation  in  A  with  analytic  coefficients  /?;(f),  satisfying 

/?j(0)  =  0,;  =  l,...,n.  (2) 

We  may  therefore  write 

0jU)  =  0'j'\  +  l3f\^  +  ....  (3) 

We  shall  restrict  f  to  a  nontrivial  real  interval  [0,eo]- 


It  is  well  known,  e.g.  [3,14],  that  the  roots  of  (1)  are  described  by  series  in  fractional 
powers  of  £.  These  series  are  commonly  called  Puiseux  series,  since  it  was  Puiseux  [23] 
who  established  their  convergence;  however  they  were  derived  formally  by  Newton  two 
centuries  earlier.  Newton,  of  course,  was  concerned  only  with  real  coefficients  and  real 
roots,  but  we  shall  consider  the  general  complex  case  first  for  simplicity.  We  shall 
obtain  the  results  we  need  by  making  use  of  a  diagram  devised  by  Newton  for  the 
purpose  of  calculating  coefficients  of  Puiseux  series  [20,24].  Since  so  many  calculation 
techniques  go  by  the  name  of  Newton,  we  shall  follow  [3]  in  calling  this  diagram  the 
Puiseux-Newton  diagram. 

Let  /3j  =  /?j'^^  be  the  first  nonzero  value  in  the  sequence  {Pj  ' ,  0j  ,■■■)■  By  def- 
inition, £j  >  l,j  =  l,...,n.  U  0j{()  is  identically  zero,  take  £j  =  oo;  also,  since  the 
coefficient  of  A„  in  P(A,  f )  is  one,  take  Co  =  0,  Po  =  I.  Now  plot  the  values  ij  versus  j, 
and  consider  the  lower  boundary  of  the  convex  hull  of  the  points  plotted.  Let  Sj  be  the 
slope  of  the  line  on  [j,j  +  1]  forming  part  of  this  boundary,  j  =  0, .  .  . ,  n  -  1.  Clearly 
1/n  <  So  <  *i  <  •  •  ■  <  «n-i-  Figure  1  shows  the  diagram  for  the  following  example 
(taken  from  [3]): 


n  =  3;Ao  =  0;  /?i(f)  =  f;  /?2(0  = -f  -  f';  Psi^) 
We  have  (o  =  0,  £i  =  1,  (7  =  I,  (3  =  2,  and  so  sq  =  «i  =  i,  • 


e^  +  2e= 


(4) 


Now  consider  the  following  "Ansatz"  argument.    Suppose  a  root  of  (1)  b  to  have 
the  form 

X(c)-Xo  =  Qc''  +  ---  (5) 

where  a  is  nonzero  and  p  is  the  smallest  power  of  c  in  the  expansion  for  this  root. 
Substituting  (5)  into  (1),  we  need 

(a"f"''+--)     +     (Ae'-  +-)(a"-'e'"-'''' +  ••■)   +   •■■ 

+     (/Jn-ie'—  +  •  •  ■)(ae''  +  •■■)  +  (/^nf'"  +  •  ■  )  =  0. 


The  terms  involving  the  smallest  powers  of  e  are  among  the  terms 

a"f"^  AQ"-'(''+'"-')P,...,/3„_iae'"-'+P,  /?„('".  (6) 

For  cancellation  to  take  place,  at  least  two  terms  with  the  same  smallest  power  of  e  must 
appear.  Equivalently,  p  must  equal  one  or  more  of  the  slopes  Sq,  . . . ,  a„_i  defined  by 
the  Puiseux-Newton  diagram.  The  following  discussion  will  apply  to  a  particular  choice 
of  such  p.  Define  /  and  d  hy  p  =  sj  =  . . .  =  sy+^-i .  so  that  the  line  in  the  diagram 
with  slope  p  passes  from  the  point  {f,(j)  to  the  point  {f  +  d,if+d)-  Cancellation  of  the 
coefficients  of  the  terms  with  the  smallest  powers  of  e  in  (6)  requires  a  to  be  the  root 
of  a  polynomial  equation  with  degree  d,  with  leading  term  /J/o''  and  constant  term 
/3y+d,  and  with  an  additional  intermediate  nonzero  term  for  each  point  {j,£j)  lying  on 
the  line  in  the  diagram  with  slope  p,  where  f  <  j  <  f  +  d.  Now  let  p  =  q/r,  where 
q,r  are  relatively  prime  integers.  By  definition,  p  is  an  integral  multiple  of  l/d,  so  d 
is  an  integral  multiple  of  r,  say  d  =  mr.  It  is  then  clear  from  the  diagram  that  of  the 
d  —  \  abcissa  values  j  between  /  and  f  +  d,  only  every  rth  value  is  a  candidate  for 
the  intersection  of  the  line  with  a  point  with  integer  coordinates.  Consequently  the 
polynomial  of  degree  d  in  a  reduces  to  a  polynomial  of  degree  m  in  a"",  which  we  may 
denote  by  Q(y).  The  conclusion  is  that  the  given  value  of  p  is  associated  with  d  roots 
with  an  expansion  of  the  form  (5),  with  a  taking  the  values 

7|'''u;-',     /j=  l,...,m,  j=  l,...r  (7) 

where  the  7^  are  the  m  roots  of  Q{j)  =  0,  7^  is  the  principal  rth  root  of  7/,  and  u>  is 
the  principal  rth  root  of  unity. 

Completing  the  example  given  above,  we  see  that  the  two  values  for  p  are  sq  —  sy  = 
^  and  S2  =  1.  In  the  case  p  =  ^  we  have  f  =  0,d  =  2,r  =  2,m  —  1,  with  Q{j)  =  7  —  1, 
so  the  possible  values  for  q  are  ±1,  giving  the  Puiseux  series 

X{()-Xo  =  ±('  +■■■■ 

In  the  case  p  =  1  we  have  f  =  2,d  =  l,r=  l,m=  1,  with  Q{j)  =  7  —  1,  so  the  only 
possible  value  for  q  is  1,  giving  the  Puiseux  series 

A(£)-  Ao  =€+■■■. 

The  subsequent  terms  in  the  series  may  also  be  calculated  by  repeating  the  process. 

The  Puiseux-Newton  diagram  may  be  used  to  establish  many  results.  For  example, 
it  is  a  trivial  consequence  of  the  diagram  that  there  is  a  Puiseux  series  (5)  with  p  =  1/n 
if  and  only  if  0n     is  nonzero. 

We  may  now  state  the  main  result  of  this  section. 

Theorem  1  Consider  the  polynomial  equation  (1),  with  roots  given  by  one  or  more 
Puiseux  series  of  the  form  (5).  Suppose  that  there  exists  Cq  >  0  such  that  all  the  roots 
A(e)  of  (1)  satisfy 

Re(A(0-Ao)<0  (8) 


Re  /?*^'  >  0, 

Re/?^'^  >0,  Im/?<'' 

=  0, 

/?(!'  =  0,    ;=3,.. 

. ,  n. 

for  0  <  f  <  eo-   Then 

(9) 

(10) 
(11) 

Proof.  The  coefficient  /?i(f)  is  the  sum  of  the  differences  Aq  -  A(f)  over  the  roots 
A(f)  of  (1);  thus  (9)  follows  from  (8),  letting  e  —  0.  The  other  results  follow  from 
the  Puiseux-Newton  diagram  as  follows.  Consider  the  Puiseux  series  corresponding  to 
p  =  So,  the  smallest  possible  value.  In  order  for  the  roots  on  the  right-hand  side  of  (5) 
to  all  be  in  the  left  half-plane,  it  is  necessary,  because  of  (7),  that  r  =  1  or  r  =  2,  with 
Re  7/i  <  0,  /i  =  1, . . . ,  m  in  the  former  case  and 

Re7h  <  0,  Im7A  =0,    h  =  \,...,m  (12) 

in  the  latter  case.  In  both  cases  p>\-  Since  the  points  ij,tj)  lie  on  or  above  the  line 
through  the  origin  with  slope  p  =  sq,  (11)  holds.  In  the  case  r  =  1,  p  >  1  so  /Jj  is 
also  zero.  In  the  case  r  =  2,  we  have 

Q(7)  =  7-+/?^iV-i+...  +  /?<'"\  (13) 


Since  -/J^^'  is  the  sum  of  the  roots  of  (5(7),  (10)  follows  from  (12).  D 

The  result  of  Theorem  1  was  first  obtained  by  Levantovskii,  using  a  completely 
different  technique  developed  for  C^  polynomial  coefficients.  We  were  unable  to  verify 
some  of  the  steps  in  his  proof,  particularly  the  justification  for  the  division  by  the 
quantity  in  [16,  p. 20,  line  D].  Although  we  do  not  doubt  that  the  result  is  correct 
for  C^  perturbations,  we  think  that  our  approach  using  the  Puiseux-Newton  diagram, 
which  requires  analyticity,  gives  far  more  insight. 

The  next  theorem  shows  that  the  conditions  (9)-(ll),  as  well  as  being  necessary  for 
(8),  are  sufficient  for  the  existence  of  a  polynomial  whose  roots  have  real  part  strictly 
less  than  Re  Aq. 

Theorem  2   Suppose  thai  (9)-(ll)  hold.    Then  there  exists  a  polynomial  of  the  form 
(I)  such  that 

Re(A(f)-Ao)  <0 

on  an  open  interval  (0,€o). 

Proof.     The  proof  is  trivial  if  n  =  1.  Otherwise,  let  l3k(i),  k  =  1, . . .  ,n,  be  defined  by 
the  coefficients  of 

(A_Ao+^f  +  r)"--(A-Ao  +  i(/?2''oU^f  +  r)(A-Ao-f(/?^''OU^e  +  e-). 
n  n  n 

(Herej  =  y^.)     D  . 


Note  that  two  roots  split  from  the  rest,  but  that  all  the  roots  have  the  same  real  part. 
(It  is  interesting  to  compare  this  with  the  more  complicated  example  in  [16,  p. 20,  line 
F],  as  well  as  with  [22,  p.  489].) 

It  is  also  possible  to  use  the  Puiseux-Newton  diagram  to  make  statements  about 
higher  order  coefficients  of  0j{c).  For  example,  it  seems  reasonable  to  suppose  that  the 
following  theorem  was  known  to  Newton,  since  its  proof  using  the  diagram  is  trivial. 

Theorem  3  Let  g  be  the  minimum,  over  all  the  Puiseux  series  describing  the  roots, 
of  the  associated  values  p.   Then 


,(t) 


0,    for  all  j  >  lk/g\. 


(14) 


Here  [x\  is  the  largest  integer  less  than  or  equal  to  x.  This  theorem  may  be  viewed  as 
a  form  of  converse  to  the  obvious  fact  that,  if /?,(f)  is  identically  zero,  j  =  n-  k, . .  .  ,n, 
then  for  any  of  the  Puiseux  series  describing  the  roots,  p  >  l/k,  since  dividing  P(A,t) 
by  the  factor  A"~*  leaves  a  polynomial  of  degree  it. 

Finally  we  note  that  Theorem  1  can  be  generalized  in  a  straightforward  way  to 
products  of  polynomials  of  the  form  (1).  Since  our  main  interest  is  in  matrices,  we 
omit  the  details  (but  see  [4]). 

3      Eigenvalues  of  Complex  Matrices 

Let  ^<°'  €  C"^"  be  a  matrix  with  eigenvalues  Ai, . . .,  A<,,  with  (algebraic)  multiplicities 
ti,...,to  respectively,  where  Re  Ai  =  . . .  =  Re  Ap  =  0  and  Re  A^+j  <  0,  .  . . ,  Re  A^  <  0. 
Let  the  Jordan  form  of  ^*°^  be  given  by 


^'0)  =  PJP-1 


where 


J  = 


Jx 


Jk 


and  the  Jordan  block 


Jk, 


A.      1 


Ja 


Jkm, 


Jkl  = 


has  dimension  riki-  We  have 

n/ti  + \-nkmi.=tk,      k-l,...,(T. 

If  mi  =  1,  A  J:  is  said  to  be  a  nonderogatory  eigenvalue,  while  if  m^  =  tk,  i.e. 
nici  =  ■  ■  ■  —  rikmi.  =  1,  ^k  is  said  to  be  nondefective  (semisimple). 

Definition.  Define  the  jth  generalized  trace  of  a  square  matrix  A,  denoted  by 

as  the  sum  of  the  elements  on  that  diagonal  of  A  which  is  j  -  1  positions  below  the 
main  diagonal.  Thus  one  obtains  the  ordinary  trace  in  the  case  j  =  1  and  the  bottom 
left  element  of  the  matrix  in  the  case  that  j  is  the  dimension  of  the  matrix.  If;  exceeds 
the  dimension  of  A,  take  tr^-'M  =  0. 

Theorem  4  Lei  A(()  be  an  analytic  matrix  function  of  c  satisfying  A{0)  =  A''°\  (By 
an  analytic  matrix  function  we  mean  a  matrix  whose  elements  are  analytic  functions.) 
We  shall  restrict  e  to  a  nontnvial  real  interval  [O.fo]-  Define  yl'^'  =  ^'(0),  the  deriva- 
tive of  A{e)  at  €  =  0.  Partition  P'M'^'P  conformally  with  the  partition  of  J  and 
denote  its  diagonal  block  corresponding  to  Jk  by  Bk,  k  =  I, .  . .  ,(t,  with  each  Bk  having 
diagonal  blocks  Bki  corresponding  <o  Jt/,  /  =  1, .  .  . ,  m^.   Define 

ri^'  =  ^tr<J'5u,    j  =  l,...,max{nt,},     k=l,...,(7. 
(=1 

(This  sum  includes  zero  terms  when  j  exceeds  riki-)  Then,  for  k  -  \, . . .  ,cr,  the  eigen- 
values of  A{()  corresponding  to  Xk  are  the  roots  of 

(Xk  -  Ao)"  +  Pki((){h  -  Ao)''-'  +  ■  ■  ■  +  Pktd^)  =  0  (15) 

where  the  0kj(()  ere  analytic,  with 

/?i;(0)  =  -^t",    J  =  l,...,max{nH). 

Proof.  Define  J{€)  =  P-^A{()P.  From  Arnold  [1,  Section  I  together  with  Theorem 
4.4],  J(e)  has  a  representation  in  terms  of  a  versal  deformation  C(e).  Specifically,  there 
exists  an  analytic  matrix  function  Y((),  with  y(0)  =  /,  such  that 

J(c)  =  Y(c)C{c)Y(c)-' 

C(c)  =  Diag(Ct(0), 

Ck(e)  =  Jk  +  Ek(e), 
and  Ek{(),  which  satisfies  Ek{0)  =   0,  is  a  "normal  form"   whose  structure  is  best 
described  by  an  example  as  follows; 


X  X     X     X  X     X 

X 

X  X  X     X 

X  X  X 

Here  we  suppose  that  Jk  has  three  blocks  JkiyJk2,Jk3,  with  dimensions  Uki  =  3, 
71^2  =  2,  nk3  =  1  respectively,  so  that  tk  =  6.  Each  X  indicates  a  complex  analytic 
function  of  e.  Denote  the  diagonal  blocks  of  Cjt(f)  by  Cki{()',  each  of  these  has  the  form 
commonly  known  as  a  companion  matrix.  Let  us  write 

1 


Ck,(e) 


'(0 


At  1 


(16) 


so  that,  for  the  example  given  above, 


Ct(0 


At  1 

Ai 

0(3).    V  q(2}. 


1 


o(^)t 


917(0     -OVi(^)     ^k-0[\'{e)  X  X 

X  Xk  1 

A' 


X 


-^IfcVCf)     A, -^-'(e)  A- 


X 


Xk-e['^(, 


Here  each  A'  indicates  a  distinct  unnamed  analytic  function  off.  Note  C/t(0)  =  Jk- 
Now  we  have 

7'(0)  =  c'(0)  +  y'(0)7  +  ,/(y-M'(0) 

=  C'(0)  +  [y'(0),7],  (17) 

where  the  latter  term  is  the  Lie  bracket 

[y"(0),j]  =  y'(0)7- jy'(o). 

Following  Arnold,  we  note  that 

([y'(o),j],/v)  =  o  (18) 

for  all  A'  such  that  A'*  commutes  with  J,  where 

{A,B)  =  tT  A'B, 


the  Frobenius  matrix  inner  product,  and  A'  denotes  the  complex  conjugate  transpose 
of  A.  Such  matrices  A''  are  block  diagonal,  K'  =  Diag(A7),  where  A'^  commutes  with 
Jk-  From  [1,9],  the  matrices  commuting  with  Jk  have  a  block  upper  triangular  Toeplitz 
form,  specifically,  for  the  example  given  above, 

b      c     d      e      f 


Kt  = 


Consequently,  A't  has  a  block  lower  triangular  Toeplitz  structure.  Note  that  if  nit  =  1 
(the  nonderogatory  case),  A't  is  a  lower  triangular  Toeplitz  matrix,  while  if  rrik  =  tk 
(the  nondefective  case),  there  is  no  restriction  on  the  structure  of  Kk-  By  (18),  then, 
we  have 

tr^-'^Git;  =  0,    j  =  l,...,max{njt,},    k=l,...,cr, 

where  the  Gki  are  the  diagonal  blocks  of  [V''(0),  J],  corresponding  in  dimension  to  Bki 
and  Jki-  (Indeed,  more  can  be  said  about  the  off-diagonal  blocks  of  [Y'(0),J]  but  we 
shall  not  need  this.)  Now,  by  definition,  the  diagonal  blocks  of  J'(0)  are  Bki,  so,  by 
(17), 

(19) 


.(^■)  =  ^tr(^)ci,(0)  =  -^^iV'(0); 


here  ^^,    is  taken  to  be  zero  if  j  exceeds  njt/. 

Since  A(()  is  similar  to  J(()  and  therefore  to  C(€),  we  need  only  examine  the 
eigenvalues  of  C(f).  We  need  the  following: 

Lemma  1 

det  (XI  -  Okie))  =  (A  -  At)''  +  /?H(f)(A  -  A^)"-'  +  •  ■  ■  +  /?,tu(0  (20) 


whe 


/?i;(0)  =  li^i'''(0)- 


In  the  case  that  At  is  nonderogatory,  the  proof  of  Lemma  1  is  trivial.  In  general,  it 
may  be  proved  by  a  fairly  straightforward  but  complicated  induction  argument  which 
gives  little  insight;  we  shall  not  give  it  here.  The  reader  may  prefer  to  verify  Lemma  1 
by  using  a  symbolic  algebra  system  such  as  Mathematica  [26]  (indeed,  this  is  how  we 
discovered  the  result).  The  proof  of  Theorem  4  is  completed  by  combining  (19)  and 
(20).  n 

The  next  theorem  continues  to  use  the  notation  of  Theorem  4. 


10 


Theorem  5   Suppose  now  thai  all  the  eigenvalues  of  A(t)  have  nonposiitve  real  part 
on  a  nontrtvial  interval  [0,io].    Then,  for  k  =  \, ...,  p, 

Rer^"<0,  (21) 

Re  r^^'  <  0,     Im  rf '  =  0,  (22) 

rt''  =  0,    >  =  3,...,max{riH}.  (23) 


and 


Proof.  First  note  that  there  is  no  restriction  on  the  eigenvalues  of  C't(f)  for  k  >  p 
as  long  as  to  is  small  enough,  since  Re  Ajt  <  0  for  k  >  p.  The  proof  is  completed  by 
applying  Theorem  1  to  equation  (15),  k  =  I, .  .  .,p.  D 

li Xk  is  nonderogatory,  i.e.  m^  =  1  and  Jk  =  Jki,  Bk  =  Bki,  we  have  r^-''  =  tr*-''Bi, 
j  =  1, . .  .  ,tk-  The  next  theorem  shows  that  the  necessary  condition  of  Theorem  5  is 
also  a  sufficient  condition  for  the  existence  of  a  stable  family  A(()  when  the  eigenvalues 
with  real  part  equal  to  zero  are  all  nonderogatory. 

Theorem  6  Assume  that  A^.  is  nonderogatory,  i.e.  mt  =  1,  for  k  =  1,  .  .  .,p.  Lei  B 
be  a  given  matrix  and  denote  its  diagonal  blocks,  partitioned  conformally  with  J ,  by 
Bie,k  =  1, . . . ,  (T.  Suppose  that  B  satisfies 

Re  tr'^'Bi  <  0, 

Re  tr '-'Bit  <  0,    Im  tr '^^Bfc  =  0, 

tr(J'5t  =  0,    j  =  3,...,«t, 

for  k  =  1, .  .  .  ,p.  Then  there  exists  an  analytic  family  A{e),  with  A{0)  =  .4'°',  A'{0)  = 
PBP~^,  such  that  the  maximum  real  pari  of  the  eigenvalues  of  A(e)  is  strictly  negative 
on  an  open  interval  {0,(o). 

Proof.  There  exists  an  Arnold  normal  form  for  J  +  eB,  namely 

J  +  eB  =  Y{e)C(i)Y(e)-\ 

where  Y{()  is  analytic  and  satisfies  Y{0)  =  I,  and  where  C(()  =  Diag(Ct(f )).  We  need 
consider  only  k  <  p  as  long  as  f  is  small  enough.  By  the  nonderogatory  assumption, 
Ci(e)  consists  of  a  single  companion  matrix  block  (see  (16)).  We  have 

B  =  C'{0)  +  [Y'{0),J]. 

Observe  that  the  relevant  diagonal  blocks  of  the  second  term  have  generalized  traces 
equal  to  zero,  so  the  first-order  terms  in  the  last  row  of  Ck{()  are  tr^^^Bk,  j  =  tk,-  ■  .,1- 
Now  define  C{e)  to  be  another  block  diagonal  matrix  with  the  same  structure  as  C{c). 
Specifically,  choose  C(()  so  that  C(0)  =  C(0)  =  J  and  C"(0)  =  C'(0),  with  the  ^-th 
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block  o{C{()  set  equal  to  the  companion  matrix  for  the  corresponding  stable  polynomial 
(of  degree  tk)  defined  by  Theorem  2.  Let 

J(e)  =  J  +  cB  +  Y{(){C{c)  -  C{())Y{c)-' 

=  Y(e)C(e)Y(()-' 

so  that  J'(0)  =  B  and,  by  Theorem  2,  J(e)  is  stable  on  an  interval  (O.to)-  The  proof 
is  completed  by  defining 

A{€)  =  PJ{()P-\     a 

We  note  that,  in  the  nonderogatory  case,  Theorem  4  is  implicit  in  the  work  of 
Levantovskii.  However,  he  apparently  did  not  obtain  our  necessary  condition  in  the 
derogatory  case.  Regarding  analyticity,  Arnold  [1,  p. 37]  mentions  that  his  normal 
form  may  be  extended  to  C^  perturbations,  although  the  proof  is  tedious.  As  in  the 
polynomial  case,  this  extension  does  not  seem  to  be  worth  the  effort  involved. 

Having  discussed  the  nonderogatory  case,  consider  now  the  other  extreme  case, 
namely  when  all  the  eigenvalues  with  real  part  equal  to  zero  are  nondefective  (semisim- 
ple),  i.e.  mic  =  tk,  k  =  1, . . . ,  p.  The  necessary  condition  of  Theorem  5  reduces  to 

Retr  Bt  <  0,      k  =l,...,p 

since  all  the  Bki  are  of  dimension  one.  This  condition  is  an  obvious  necessary  condition 
for  stability,  and  equally  clearly  not  sufficient  for  the  existence  of  a  stable  family  A{e) 
with  v4'(0)  =  PBP~^ .  A  necessary  and  sufficient  condition  is,  however,  well  known  in 
this  case,  as  the  following  theorem  states. 

Theorem  7  Assume  that  Xk  ts  nondefective,  i.e.  rrik  =  tk,  for  k  =  1, . .  .,p.  Lei  A{e) 
be  any  matrix  family  with 

A(0)  =  yl<°',  A'(0)  =  PBP-\  (24) 

and  where,  as  before,  Bk  is  the  diagonal  block  of  B  corresponding  to  \k-  Then  a 
necessary  condition  for  all  the  eigenvalues  of  A(()  to  have  nonposilive  real  part  on  a 
nonirtvial  real  interval  [0,  fo]  [respectively  sufficient  condition  for  all  the  eigenvalues  of 
A{()  to  have  negative  real  part  on  an  open  interval  (0,  (q)]  is  that,  for  k  —  1, . . . ,  p,  the 
maximum  real  part  of  the  eigenvalues  of  Bk  is  less  than  or  equal  to  zero  [respectively 
less  than  zero]. 

The  proof  uses  [14,22].  Note  that  in  this  case  the  relevant  first  p  columns  of  P 
are  eigenvectors.  Note  also  that,  unlike  in  Theorem  6,  the  sufficient  condition  applies 
to  any  family  >4(f)  satisfying  (24);  this  is  because  the  eigenvalues  are  Lipschitz  in  t. 
Analyticity  of  A(()  is  not  required;  C^  parameter  dependence  is  sufficient. 

When  Ajt  is  simple,  i.e.  rrik  =  tk  =  I,  k  =  1, . . .  ,p,  the  necessary  conditions  of 
Theorems  5  and  7  and  the  sufficient  condition  of  Theorem  6  all  reduce  to  the  following: 
the  first  p  diagonal  elements  of  B  have  nonpositive  real  part. 

The  most  difficult  case  is  therefore  the  derogatory,  defective  case,  for  which  we  have 
given  only  a  necessary  condition.  Perhaps  a  sufficient  condition  could  be  derived  by 
consideration  of  the  off-diagonal  blocks  of  [y'(0),  J]. 
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4     The  Real  Case 

Every  eigenvalue  of  a  real  matrix  is  either  real  or  one  of  a  complex  conjugate  pair.  The 
results  given  in  Section  3  apply,  in  particular,  in  the  real  case.  However,  they  may  be 
simplified.  When  /1*°'  is  real,  its  Jordan  form  may  be  taken  to  have  the  form 


where  P^Q  =  I  and 


A^°^  =  PJQ^ 

J  =  Diag(Jo,  Ji,7i,...,y,,7,,...); 
P=[Po,Pi,Pi,...,P„P,,...]; 

Q  =  [Qo,Qi,Qi,..-,Q,,Q, ,■■■]■ 

Here  Jo  corresponds  to  the  eigenvalue  Xq  =  0,  the  only  possible  real  eigenvalue  on  the 
imaginary  axis;  this  block  and  the  corresponding  blocks  Pq,  Qq  may  be  absent.  The 
blocks  Pq  and  Qo  can  be  taken  to  be  real.  Blocks  J\, . . . ,  J ,  correspond  to  eigenvalues 
with  positive  imaginary  part  and  zero  real  part.  The  symbol  1,  of  course,  denotes  the 
complex  conjugate  of  ;.  The  number  of  eigenvalues  on  the  imaginary  axis  is  p  =  2s+  1 
if  the  zero  eigenvalue  is  present  and  p  =  2s  otherwise.  Blocks  not  explicitly  listed  are 
of  no  interest  since  their  eigenvalues  have  negative  real  part. 

Now  assume  also  that  A(i)  is  real.  We  have,  if  the  zero  eigenvalue  is  present,  that 
Bq  =  QIA'^'^-^Pq  is  real,  and  hence  that  Tq  is  real  for  all  j.  Consequently,  conditions 
(21)-(23)  reduce,  in  the  case  it  =  0,  to  r^^^  <  0,  r^^'  <  0,  and  t'^J^  -  0  for  j  >  2. 

Now  consider  a  complex  conjugate  pair  of  eigenvalues.  In  this  case  the  t^,  are 
generally  complex,  but  they  occur  in  complex  conjugate  pairs.  Since  each  of  condi- 
tions (21)-(23)  applies  equally  to  a  complex  number  or  its  conjugate,  each  condition 
appears  exactly  twice.  Thus  it  is  necessary  to  apply  them  only  to  Tjf'  corresponding 
to  eigenvalues  with  positive  imaginary  part. 

Most  likely  it  is  possible  to  rewrite  these  conditions  using  only  real  quantities  by 
using  the  real  Jordan  form  [8,15].  However,  our  attempts  to  do  this  have  led  only  to 
complicated  results.  See  [8]  for  the  extension  of  the  Arnold  normal  form  to  the  real 
field. 

Regarding  complex  conjugate  roots  of  polynomials  with  real  coefficients:  instead  of 
(1)  we  may  consider 

P(A,0     =     [(A-ia)"+/?i(f)(A-fa)"-^  +  ...  +  /?„(0]   • 

[(A  +  ia)"  +  0i(()(X  +  ia)"-'  +  •  ■  •  +  ^„(0]  (25) 

where  a  is  real  and  nonzero.  Suppose  now  that  P(X,()  has  no  root  with  positive  real 
part  on  a  nontrivial  interval  [0,eo]  By  the  generalization  of  Theorem  1  mentioned  at 
the  end  of  Section  2,  it  follows  that  (9)-(ll)  hold.  Let 

R{X)  =  X^  +  a'^. 
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Multiplying  out  the  factors  in  (25)  we  see  that 

P(X,  e)  =  fl(A)"-2(fl(A)2  +  uieXRiX)  +  U2€R{X)  +  use)  +  O(i^) 

where  i/i  >  0,  ^3  <  0,  and  there  is  no  sign  restriction  on  u^.  The  same  result  (with  a 
different  derivation)  is  given  in  [16,  p. 20,  line  E]. 

5      Reducing  the  Spectral  Radius 

The  spectral  radius  of  a  matrix  is  the  radius  of  the  smallest  disk  in  the  complex  plane 
containing  all  its  eigenvalues.  For  the  purposes  of  this  section,  change  the  definition  of 
matrix  stability  to  require  the  spectral  radius  to  be  less  than  one.  We  first  modify  the 
polynomial  theorems  accordingly. 

Theorem  8  Consider  the  polynomial  equation  (1),  with  roots  given  by  one  or  more 
Puiseux  series  of  the  form  (5).  Suppose  that  there  exists  €0  >  0  such  that  all  the  roots 
A(f )  of  (1)  satisfy 

\X(c)\  <  |Ao|  (26) 

for  0  <  f  <  fo-   Then 

ReAo/?l''>0,  (27) 

Re  Aq/?^^'  >  0,  Im  Ao4^'  =  0,  (28) 

/?]''  =  0,    ;=3,...,n.  (29) 

Proof.  First  note  that  when  Ao  is  real  and  positive  the  result  is  the  same  as  before; 
this  is  because  the  circle  of  radius  |Ao|  centered  at  the  origin  lies  to  the  left  of  the 
vertical  axis  Re  r  =  Re  Aq-  When  Aq  is  real  and  negative  we  obtain  the  negative  of 
the  original  conditions,  since  the  circle  lies  to  the  right  of  the  corresponding  axis.  In 
general,  (26)  is  equivalent  to 

ReAo(A(f)-Ao)<-^   |A(0  -  Ao|   ' 

so  that  a  necessary  condition  for  (26)  to  hold  is  that 

ReAo(A(f)-Ao)<0.  (30) 

The  proof  is  now  an  easy  modification  of  the  proof  of  Theorem  1  as  follows.  The 
coefficient  /3i(f)  is  the  sum  of  the  differences  Aq  —  A(e)  over  the  roots  A(f)  of  (1);  thus 
the  first  condition  follows  from  (26),  letting  €  —  0.  The  other  results  follow  from  the 
Puiseux-Newton  diagram.  In  order  for  the  roots  on  the  right-hand  side  of  j_5)  to  lie  in 
the  half-plane  defined  by  (30),  it  is  necessary  that  r  =  1  or  r  =  2,  with  Re  Ao7h  <  0  in 
the  former  case  and 

— 2  — 2 

Re  Ao7h  <  0,  Im  AqT^  =  0,  /i  =  1, . . . ,  m 


14 


in  the  latter  case.  The  rest  of  the  proof  is  exactly  as  before.  O 

A  sufficient  condition  is  somewhat  harder  to  obtain  in  the  spectral  radius  case. 
The  difficulty  is,  by  way  of  example,  that  an  eigenvalue  splitting  of  the  form  1  ±  jfi , 
while  having  the  same  real  part  as  the  unperturbed  eigenvalue  1,  has  spectral  radius 
\/\  +  c.  Consequently,  we  need  to  make  an  assumption  that  the  left-hand  side  of  (27) 
is  sufficiently  bigger  than  the  first  left-hand  side  of  (28).  For  example,  we  have: 

Theorem  9  Suppose  that  Aq  ?^  0  and  that  (27)-(29)  hold.  Suppose  also  that 


\p\"\  <2ReAoi^.  •  (3i; 


jIDi   ^  oD„T  ^1 


Then  there  exists  a  polynomial  of  the  form  (1)  such  that 

|A(e)|  <  |Ao| 
on  en  open  interval  (0,  to)-     < 
Proof.     Let  /^it(f),  t  =  1, . . . ,  n,  be  defined  by  the  coefficients  of 

^(1)  ^(1)  rf\) 

(A  -  Ao  +  ^f)""'(A  -  Ao  +  i(4"0^  +  ^0(A  -  Ao  -  i{p['\)^^  +  ^e).    □ 

Note  that  the  larger  |Ao|  is,  the  more  likely  (31)  is  satisfied  for  given  p\  ,  /?2  •  This 
reflects  the  fact  that  the  the  smaller  a  circle  is,  the  greater  the  distance  between  it 
and  a  tangent  line  segment  of  given  length.  In  the  limit  case  Aq  =  0,  reduction  of  the 
spectral  radius  is,  of  course,  impossible.  The  factor  n  in  (31)  reflects  the  particular 
choice  of  polynomial  used  in  Theorem  2  and  can  be  improved  (see  [4]). 

Now  take  A^^^  to  have  the  same  Jordan  form  as  in  Section  3,  but  assuming    |Ai|  = 
. . .  =  |Ap|  =  1  and    |Ap+i|  <  1, . . . ,   jA,,]  <  1.  Instead  of  Theorem  5,  we  obtain: 

Theorem  10  Suppose  that  the  spectral  radius  of  A(e)  ts  less  than  or  equal  to  one  on 
a  nontrtvtal  interval  [0,  fo]-    Then,  for  k  =  I, . . . ,  p, 

Re  Atr^^'  <  0, 

ReA^rf'<0,     ImA^rf  =  0, 

and 

r^-"  =  0,     j  =  3 maxjni;}. 

The  proof  is  identical  to  the  proof  of  Theorem  5,  using  Theorem  8  instead  of  Theorem  1. 
In  the  nonderogatory  case,  we  obtain  the  following  sufficient  condition. 
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Theorem  11  Assume  that  A*  is  nonderogaiory,  i.e.  ruk  =  1,  for  k  =  \,...,p.  Let 
B  be  a  given  mairtx  and  denote  its  diagonal  blocks,  partitioned  conformally  with  J ,  by 
Bk ,  k  =  1 , .  . . ,  <T.  Suppose  that  B  satisfies 

ReAttr<''Bi  <  0, 

Re  A^tr'^'Bi  <0,    Im  A^tr '^^B*  =  0, 
tr(^>Sjt  =0,    j  =  3,...,tk, 
for  k  =  I, . .  .,p.  Suppose  also  that 

|tr(2)Si|  < --ReAitr'^'Bfc.  (32) 

Then  there  exists  an  analytic  family  A{c),  with  A(0)  =  A^°\  A'{0)  -  PBP-\  such 
that  the  spectral  radius  of  A(c)  is  strictly  less  than  one  on  an  open  interval  (0,€o)- 

The  proof  combines  those  of  Theorem  6  and  9.  As  before,  (32)  can  be  weakened. 

Finally,  in  the  nondefective  case  we  have  the  following  modification  to  Theorem  7. 

Theorem  12  A  necessary  condition  for  the  spectral  radius  of  A(e)  to  be  less  than  or 
equal  to  one  on  a  nontrivial  real  interval  [0,  to]  [respectively  sufficient  condition  for 
the  spectral  radius  of  A(()  to  be  less  than  one  on  an  open  interval  (0,eo)7  ts  that,  for 
k=  l,...,p,  the  maximum  real  part  of  the  eigenvalues  of  XkBk  is  less  than  or  equal  to 
zero  [respectively  less  than  zero]. 

The  proof  uses  [22]. 

6      Practical  Implications 

How  useful  are  our  results  in  practice?  The  necessary  condition  is  practical  in  the  sense 
that,  if  the  Jordan  form  of  ^4'°'  is  known,  a  candidate  first-order  perturbation  A^^^  can 
be  rejected  immediately  if  it  does  not  satisfy  the  given  necessary  condition.  However, 
the  sufficient  condition  is  of  little  practical  use,  since  nothing  is  said  about  how  to 
compute  the  actual  matrix  family  A{e).  Indeed,  it  is  well  known  that  the  nonsymmetric 
eigenvalue  problem  is  very  ill-conditioned  in  the  case  of  multiple  eigenvalues.  The 
reason  for  the  ill-conditioning  is  immediately  apparent  from  the  Puiseux  series  (5);  the 
eigenvalues  are  not  Lipschitz  when  the  Puiseux  exponent  is  less  than  one.  Computing 
the  Jordan  form  of  yl'°'  is  itself  a  very  ill-conditioned  numerical  problem  and,  in  fact, 
it  is  the  instability  of  this  process  which  originally  led  to  Arnold's  normal  form  for 
perturbed  matrices. 

Several  proposals  have  been  made  to  stabilize  the  computation  of  the  Jordan  form 
of  a  matrix  ^'°';  see  [6,10,13].  The  point  has  often  been  made  that  a  block  Schur  form 
(a  block  diagonal  form  with  a  full  upper  triangular  block  for  each  eigenvalue  [11,  p. 338]) 
is  much  more  stable.   The  idea  of  the  Schur  form,  eis  opposed  to  the  Jordan  form,  is 
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to  avoid  the  reduction  to  too  small  a  number  of  parameters.  Unfortunately  this  is  of 
little  use  for  our  purposes  since  the  Arnold  theory  requires  versa!  deformations,  which 
involve  using  a  minimal  number  of  parameters;  hence  the  use  of  the  Jordan  form  and 
the  block  triangular  Toeplitz  matrices  which  commute  with  the  Jordan  form.  Somehow 
one  wishes  to  compute  the  "least  generic"  of  the  Jordan  forms  of  those  matrices  which 
are  within  some  tolerance  of  A^°\  Perhaps  the  most  hopeful  approach  along  these 
lines  is  that  of  Fairgrieve  [7],  who  has  introduced  a  path  following  method  for  this 
purpose.  Fairgrieve's  method  is  based  on  Arnold's  normal  form  and  in  particular  uses 
Arnold's  classification  of  the  degeneracy  hierarchy  [1],  in  which,  as  is  well  known,  a 
nonderogatory  eigenvalue  is  the  most  generic  and  a  nondefective  one  the  least  generic, 
but  for  which  the  ordering  of  the  "middle"  cases  is  by  no  means  obvious. 

We  conclude  with  the  description  of  two  numerical  experiments,  conducted  in 
Matlab  [19],  which  give  substantial  insight.  In  the  first  experiment,  we  set  J  to  the 
Jordan  block  of  order  n  =  10  with  zeros  on  the  diagonal  and  ones  on  the  superdiagonal. 
We  then  computed  the  spectral  radius  of  J  +  fB<*',  where  e  =  10~^°  and  B^*^  is  a 
random  real  matrix  of  order  n  with  the  property  that 

tr'-'^S"^' =  0,     j  =  n-k+l,...,n.  (33) 

Thus,  for  k  =  1,  the  bottom  left  element  of  B<*'  is  zero;  in  addition,  for  k  =  2,  the  two 
elements  on  the  (n-  l)th  subdiagonal  sum  to  zero,  etc.  The  results  were  averaged  over 
100  random  matrices  and  the  resulting  averaged  spectral  radii  plotted,  as  a  function 
of  k,  as  "*"'s  in  Figure  2.  The  "+"  symbols  in  Figure  2  show  the  quantities  fi/("-*). 
Inspecting  Figure  2,  one  sees  a  steady  reduction  in  the  averaged  spectral  radius  as 
k  increases  from  0  to  n/2,  with  close  agreement  between  the  two  plotted  values.  At 
k  =  n/2,  however,  a  "floor"  is  reached.  Variation  off  and  the  machine  precision  showed 
that  this  floor  is  not  a  consequence  of  round-off  error,  as  might  be  suspected.  Instead, 
it  is  easily  explained  by  the  Puiseux-Newton  diagram.  By  imposing  (33),  we  are  setting 
the  corresponding  polynomial  coefficients  P,  to  zero.  Thus  as  k  increases  from  0  to 
n/2,  slopes  with  value  1/n,  l/(n  -  1), . . . ,  l/(n/2  +  1)  are  eliminated  as  possibilities. 
Consequently,  roots  of  the  form  (5)  with  corresponding  powers  are  not  possible,  and  we 
see  the  close  association  of  the  "*"  and  "+"  plots.  However,  when  k  reaches  7i/2,  the 
second-order  coefficient  /?„  comes  into  play.  Since  the  second-order  terms  arising  from 
the  Arnold  normal  form  are  not  zero,  the  slope  with  value  2/n  cannot  be  eliminated 
from  the  Puiseux-Newton  diagram.  Hence  the  "floor".  We  have  succeeded  only  in 
reducing  the  order  of  the  eigenvalue  perturbations  from  e^'"  to  e^^". 

In  the  second  experiment,  we  set  7  to  a  full  Jordan  block  of  order  n  and  defined  B 
to  be  a  random  matrix  with  the  property  that 

tr'^'S  =  0,     j=  1;    j  =  3,...,n. 

This  time  we  plotted  the  maximum  real  part  of  the  eigenvalues  of  J  +  tB  against  the 
value  of  tr  '^'5.  Again,  the  results  were  averaged  over  100  random  cases.  Plots  are 
shown  in  Figure  3  for  n  =  4  and  n  =  5.  The  reason  for  the  two  very  different  graphs 
is  similar  to  the  explanation  for  the  result  of  the  first  experiment.   In  the  case  n  —  b, 
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the  nonzero  second-order  coefficient  /Jj  '  yields  a  slope  with  value  2/5  in  the  Puiseux- 
Newton  diagram;  this  prevents  a  reduction  in  the  maiximum  real  part  of  the  perturbed 
eigenvalues.  However,  in  the  case  n  =  4,  the  slopes  in  the  diagram  must  be  at  least 
1/2,  so  that  stability  is  possible  when  tr  ^^'B  <  0.  Theorem  5  suggests  that  the  larger 
the  quantity  — tr  *-'B  is,  the  "more  stable"  the  perturbed  matrix  should  be,  and  indeed 
this  is  supported  by  the  graph  for  n  =  4.  Stability  is  not  actually  achieved,  however, 
because  of  the  effect  of  f3^'    which  is,  in  general,  nonzero. 
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